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ABSTRACT 

Semiconvection - mixing that occurs in regions that are stable when considering compositional gra¬ 
dients, but unstable when ignoring them - is shown to have the greatest potential impact on main 
sequence stars with masses in the range 1.2 Mq — 1.7 Mq. We present the first stellar evolution 
calculations using a prescription for semiconvection derived from extrapolation of direct numerical 
simulations of double-diffusive mixing down to stellar parameters. The dominant mode of semicon¬ 
vection in stars is layered semiconvection, where the layer height is an adjustable parameter analogous 
to the mixing length in convection. The rate of mixing across the semiconvective region is sensitively 
dependent on the layer height. We find that there is a critical layer height that separates weak semi¬ 
convective mixing (where evolution is well-approximated by using the Ledoux criterion) from strong 
semiconvective mixing (where evolution is well-approximated by using the Schwarzschild criterion). 
This critical layer height is much smaller than the minimum layer height expected from simulations 
so we predict that for realistic layer heights, the evolution is nearly the same as a model ran with the 
Schwarzschild criterion. We also investigate the effects of compositional gradient smoothing, finding 
that it causes convective cores to artificially shrink in the absence of additional mixing beyond the 
convective boundary. Layered semiconvection with realistic layer heights provides enough such mix¬ 
ing that stars will still evolve as if the Schwarzschild criterion is employed. Einally, we discuss the 
potential of detecting such semiconvection and its implication on convective core sizes in solar-like 
oscillators. 

Subject headings: stars: evolution — stars: interiors — convection 


1. INTRODUCTION 

One of the outstanding challenges in constructing accu¬ 
rate stellar models is understanding macroscopic mixing 
driven by fluid instabilities. Direct numerical simula¬ 
tions of fluid instabilities, however, cannot be included 
in stellar evolution calculations over an appreciable frac¬ 
tion of a star’s life, due to the large separation of length 
and time scales between small-scale fluid motion and 
global stellar behavior. Instead, stellar evolution codes 
rely on one-dimensional prescriptions for the transport of 
energy, composition, and angular momentum associated 
with instabilities such as convecti on, overshoot, shear , 
etc. Mixing length theory (MLT) (Bohm-Vitense 19581 
is the standard way of modeling convection m stars and 
planets, and is usually applied in two steps. First, one 
determines which regions are convectively stable and un¬ 
stable, then one adds a model for the convective flux to 
the total heat flux as well as a model for turbulent mixing 
in the evolution equation for each chemical species. The 
first of these two steps is typically done in a local man¬ 
ner through a comp arison of structural gradie nts within 
the object (see also |Kippenhahn et al.||2012|). Ignoring 
composition gradients, the presence or absence of convec¬ 
tion is established by the Schwarzschild criterion, which 
states that a region is stable when 
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where Vrad is the radiative temperature gradient (the 
temperature gradient required if all the energy is trans¬ 
ported via photon diffusion), Vad is the adiabatic tem¬ 
perature gradient, computed directly from the equation 
of state, K is the opacity, P is the pressure, I is the lumi¬ 
nosity, a is the radiative constant, c the speed of light, G 
the gravitational constant, T the temperature, and m the 
local mass coordinate. To account for composition gradi¬ 
ents, one must use the Ledoux criterion instead, whereby 
stability occurs when 


Vrad < Vad + ^V^ = V^, (2) 

where (j) and 5 are thermodynamic derivatives of the 
equation of state. 


and 
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is the local non-dimensional composition gradient, with fj, 
being the mean molecular weight of the material. We will 
refer to stellar models computed using the Schwarzschild 
criterion as the Schwarzschild case or Schwarzschild mod¬ 
els, and models computed using the Ledoux criterion will 
be similarly named. 

If the composition gradient is stabilizing (V^ > 

0, so Vad < Vl) then there can be regions where 
Vad < Vrad < Vl that are stable to convection under the 
Ledoux criterion, but unst able under the Schwarzschild 
criterion. It is now known (Kato|1966 1 that such regions 
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can be subject to a different kind of instability called os¬ 
cillatory double-diffusive convecti on (ODDC) and are of¬ 
ten referred to as semiconvective (Schwarzschild & Harm 


19581. Being diffusive in nature, semiconvective mixing is 
thought to be much weaker than convective mixi ng, and 
is typically modeled with a separa te prescription (Langer^ 
et al.|19^ Castellani et al.|1985 ), the strength of which 
is controlled by additional paraineters analogous to the 
mixing length in convection. 

Semiconvection was his torically used in modeling high- 
mass stars > 10 MfT) (ISchwarzschild &: Harm 1958 


Stothers 1970[ Stothers &: (Jhin|1975[ hanger et al. 




Semiconvective mixing has been in vestigated tor lower- 

mass stars on the main sequence ([Faulkner fc Cannon 


19731 Gabriel fc Noels|1977 [Crowe Mitalas|1982 ' 5iTTO 


Aguirre et al. 2010b|), but not using prescriptions moti- 
vated by numerical simulations. Semiconvective mixing 
in main sequence stars is typically encountered during 
convective core burning, just outside of the region that 
is Ledoux-unstable to convection. It is caused by the 
buildup of composition gradients induced by the long 
tail of low-temperature proton-proton chain burning that 
extends outside of the convective core for main sequence 
star s in the mass range 1.2 Mq — 1.7 Mq (see section 
4.11. Generally speaking, they occur whenever there is a 
stabilizing composition gradient adjacent to a convection 
zone, such as near H-burning shells in evolved stars as 
well as in later stages of massive star evolution. 

In this paper we implement the new prescript i on for 
semiconvective mixing proposed by Wood et al. (|2013 l 
into the open-source stellar evolution code MESA^ Pax- 


ton et al. 2011 2013 ). We review the general properties 
ot the Wood et al.j model in section We argue that 
ODBC most likely takes the form of layered semiconvec¬ 
tion fomparameters representative of stellar interiors (see 
also Garaud|2014[ ) , with the layer height remaining a free 
parameter. Sectionj^details how we implement this mix¬ 
ing prescription in MESA. It also discusses the effects of 
numerical smoothing on convective core evolution. Sec¬ 
tion presents applications of this theory to the growth 
of convective cores in main sequence stars and identi- 
hes the stellar mass range which is most sensitive to the 
effects of semiconvection. We show the effects of vary¬ 
ing layer heights on convective core growth, and derive a 
critical layer height for semiconvective regions to persist 
throughout the main sequence. We summarize our con¬ 
clusions in section]^ and briefly outline future directions 
of investigation into semiconvective mixing. 

2. OSCILLATORY DOUBLE-DIFFUSIVE 
CONVECTION 

As discussed in section oscillatory double diffusive 
convection is a mild form of convection that occurs in the 
presence of stabilizing chemical gradients and destabiliz¬ 
ing thermal gradients. The onset and strength of ODDC 
is controlled by several dimensionless parameters. The 
first is the Prandtl number, 


Pr = —, 
f^T 


(5) 


which is the ratio of the microscopic kinematic viscosity 
V to the thermal diffusivity kt (both with units of cm^/s 
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in cgs units). Low Prandtl numbers (« 10“® — 10“^) 
are typical in stars and gas giant planets, while telluri c 
planet interiors may have Pr > 1 ( Soderlund et al.|2013 ). 
The second parameter is the diffusivity ratio. 


Hi'j' 


( 6 ) 


of the microscopic compositional diffusivity to the 
thermal diffusivity. ODDC only occurs when r < 1, 
which is the standard situation in stars since heat is 
transported via photon diffusion as well as collisions be¬ 
tween nuclei, while chemical species only diffuse through 
collisional processes. Finally, the third parameter is the 
density ratio 

~ -’ O) 

where 6 and (j) are the equation of state derivatives de- 
hned in the previous s ection (with S = (j) = 1 for an ideal 
gas). As discussed by Baines & Gill (19691, semiconvec¬ 
tion occurs when 


where 


-^crit ^0 1; 

(8) 

Pr 4- r 

- Pr + 1 ■ 

(9) 


Equation ([8| is a local criterion, analogous to the 
Schwarzschild and Ledoux criteria used to determine con¬ 
vective instability. It can be rewritten in terms of V as 


Vad + Rcrif^V^ <V <Vad + '^V^. 


( 10 ) 


For Pr,T ^ 1 we have Rent ^ Ij so the semicon¬ 
vective criterion here is nearly what is typically em- 
ployed in stars, Vad < V < Vad + fV^ (Langer et al. 

19851. Regions with Rq > 1 {V > Vad + f 
unstable to convection, while regions with Rg < i?crit 
(V < Vad + RcritfVfj) are stable and therefore radia¬ 
tive. 

The region of parameter space unstable to ODDC is 
itself divided into two distinct domains, one in which the 
semiconvection later transitions into a layered state, and 
one in which it does not. The former has been deter- 


mined empirically through simulations by Mirouh et al. 
(20121, to encompass virtually all of the ODDC region 


for stellar parameter regimes of Pr, r 1. In what fol¬ 
lows, we therefore ignore the possibility of non-layered 
semiconvection. Layered semiconvection takes the form 
of fully convective layers of height 77^, separated by thin, 
stably stratihed interfaces. For stellar parameters, the in¬ 
terfaces are very mobile, and tend to merge rapidly with 
one another. As a result, the ultimate layer height in 
semiconvection is difficult to determine from numerical 
simulations and will be treated as a free parameter of the 
model, much like the mixing length in standard MLT. 

The total heat flux at a given location in a star can be 
expressed as the sum of the radiative and semiconvective 
heat fluxes, 

Ftot = T^rad + ^semi; (11) 

where the radiative heat flux is given by 


-^rad — ^rad » ? 

dr 


( 12 ) 
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and the semiconvective heat flux is given by 
= -A:rad(NuT - ^ ^ 


dr dr 


ad 


(13) 


where = pcpKT is the microscopic thermal conduc¬ 
tivity and Nut is the thermal Nusselt number (see be¬ 
low). Similarly, the turbulent compositional flux is writ¬ 
ten, 


Ff, = -k^(Nu^- 1)^, 


(14) 


where Nu„, is the compositional Nusselt number. Wood 
et al. (2013|) ran a systematic set of direct numerical sim¬ 
ulations of layered semiconvection to measure the Nusselt 
numbers Nut and Nu^ as functions of the model param¬ 
eters Ro, Pr, r, and Hp. They found that both Nusselt 
numbers vary only weakly with Rq, and depend on Hl 
only through the Rayleigh number, Ra, given by 

(15) 

(16) 
(17) 

The next section details how this mixing prescription is 
actually implemented in MESA. 
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Ra = 

Empirical fits to the data suggest that 

Nut - 1 = 0.1 Pr^/^Ra^/^, 

and 

Nu^ - 1 = 0.03 


3. IMPLEMENTATION IN MESA 

3.1. Semiconvective mixing prescription 

MESA is a one-dimensional Lagrangian code that 
solves the equations of stellar structure in a fully coupled 
(unsplit) method, although options to split the mixing 
equations from the burning and structure equations ex¬ 
ist. The two equations that are directly affected by the 
addition of extra mixing are the equations for thermal 
and compositional transport. The Hrst is discretized as 


Tfc-i — Tk 


dmk 



Tk 

\dm ) 

Tk_ 


(18) 


where Tj, is the temperature at the center of cell k, 
is the mass-interpolated temperature at the outer face 
of cell k (similar for pressures Pk and Pk), dmk is the 
mean mass of cells k and k — 1, and V/j is the tempera¬ 
ture gradient at the face of cell k (cell indices in MESA 
increase inward towards the core). While not explicit 
in this equation, the semiconvective prescription directly 
affects the the calcu lation of the local temperature gra¬ 
dient, V (see section 3.1.1 below). The equation of com¬ 
positional transport on the other hand, is 

Xi^k-i{t + St) — Xi^k{t) = —(Ti,fc+i ~ Fi^k) -— , 

dt dmk 

(19) 


where St is the time step and 

Fi^k = {Xi^k ~ Xi^k-i) 


Dk 

dmk 


( 20 ) 


is the flux of the i*** species through the outer face of the 
cell, Dk is the compositional diffusion coefficient, and 


Xi^k is the mass fraction of the sp ecies in the cell, 
see section 6.2 in Paxton et al. (2011) for a more detailed 


discussion. In this case, semiconvective mixing directly 
affects the calculation of Dk- 
The mixing type of a given cell is determined entirely 
locally, with regions labeled as semiconvective when 


^ad ^ ^rad ^ ^ad “t” 


( 21 ) 


3.1.1. Temperature gradient calculation 

As in standard MLT, we use the relationship between 
the heat flux and V to solve for the latter. From equation 
(13), we rewrite the semiconvective flux as. 




Tp^cpgn^' 


P 


P 


(V - Vad) Hi 


1/3 


(V — Vad) , 
( 22 ) 


where cp is the specific heat at constant pressure, and 
g = Gm/r"^ is the local gravitational acceleration. Di¬ 
viding equation (11) by kj-i^^Tpg/P yields, 

(23) 

(24) 


Vrad — V -|- Vs 


where 


„ _ -Fsemi f dlnPV 

V semi — 1 rr \ 1 I 

^rad^ \ dr J 


Tpgha.d 


Equation (23) is a fourth-order polynomial equation in V 


so could in principle be solved analytically for the latter. 
However, the expression of the solution is too compli¬ 
cated to be of practical value. Instead, we solve it nu¬ 
merically as part of the overall Newton solve to advance 
the star in time. The Jacobian for this Newton solve also 
requires partial derivatives of V with respect to local stel¬ 
lar variables, such as P, T, m{r), etc. We calculate these 
by differentiating equation (23) implicitly. 


3.1.2. Compositional mixing coefficient calculation 

The compositional mixing rate is locally determined by 
the semiconvective diffusion coefficient, Psemi- Its value 
in layered semiconvection is given by 


Usemi — 1) 


^0.38 

= 0.03^Inr 


Spg'^ 


1X7 . 


u4 


0.37 


(25) 


and is calculated once V is known. As with V, partial 
derivatives of Psemi are necessary for the J acob ian, which 
are computed by differentiating equation (25). 


3.2. Effects of composition gradient smoothing on 
convective boundaries 

When using the Ledoux criterion without additional 
mixing schemes such as semiconvection or overshoot, the 
resulting composition gradients are often not smooth. By 
default, MESA smooths using a Gaussian smoothing 
formula with a 7-cell window. This kind of smoothing 
can have significant, but artihcial effects on stellar evo¬ 
lution. Indeed, in core-convective stars, there is usually 
a sharp jump in composition between the outermost con¬ 
vective cell and the cell immediately above it. This pro¬ 
duces a spike in which grows in time, and is the main 
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Fig. 1. — Evolution of Vi profiles of a 1.5 Mq star when smooth¬ 
ing is suddenly turned on during evolution on the main sequence. 
Profiles of Vi as a function of mass coordinate are shown for five 
consecutive models where V^-smoothing is turned off (red), as well 
as five subsequent consecutive models with V^-smoothing enabled 
(blue). The extent of the convective core is shown by the dashed 
lines and is constant in mass coordinate until V^-smoothing is 
turned on, after which the outer mass coordinate of the convective 
core moves inward as the compositional gradient pushes into the 
formerly convective region. 

reason why in the absence of smoothing or any kind of 
mixing beyond the convective boundary, convective cores 
do not grow as large in the Ledoux case as compared to 
the Schwarzschild case. If this spike is smoothed before 
the MLT module computes the convective boundaries, 
then the region of high-V^ is artificially spread inward 
into cells that used to be convective. This then causes 
the convective boundary to move inward relative to that 
of a run which performs the same steps without compo¬ 
sitional gradient smoothing. 

Figure shows the evolution of a compositional dis¬ 
continuity at the convective boundary in a numerical ex¬ 
periment in which compositional smoothing is suddenly 
switched on during the main sequence. This clearly 
illustrates how the spike in that had built up at 
the convective boundary starts to diffuse inward when 
smoothing is turned on, causing the convective bound¬ 
ary to artificially move inward with it. Figure!^ shows a 
Kippenhahn diagram of a similar run, as weTT as that 
of a star where no smoothing was used. As soon as 
smoothing is enabled, the convective core starts mov¬ 
ing inward, resulting in a significantly smaller convective 
core for the remainder of the main sequence compared 
to the run without smoothing. This effect occurs for all 
stars that develop a convective core - not just stars that 
show a strong difference between their evolution under 
the Ledoux and Schwarzschild criteria; smoothing al¬ 
ways shrinks a convective core relative to the same model 
without smoothing. We consider this to be an artificial 
effect and our models are run without V^-smoothing, 
unless otherwise noted. We will see in section |4] that this 
distinction is rendered moot when layered semiconvec¬ 
tion is accounted for because of the mixing that naturally 
occurs beyond the convective boundary. 



Fig. 2. — Comparison of the impact of V^-smoothing on con¬ 
vective core evolution of a 1.5 Mq star. The top panel shows 
a Kippenhahn diagram of evolution under the Ledoux criterion 
without V^-smoothing, where the vertical lines mark the extent 
of convective regions at each time step. The bottom panel shows 
the evolution of the same star, but with V^-smoothing turned on 
suddenly (as in Figure Q when the core hydrogen mass fraction 
drops below 0.4 which occurs at £ = 9.5 x 10® yrs. As soon as V^^- 
smoothing is turned on, the core begins to shrink at a near constant 
number of cells per step. The apparent delay in core shrinkage in 
mass coordinate is merely due to the increased spatial resolution 
near the convective boundary. 

4. APPLICATIONS TO STELLAR EVOLUTION 

4.1. Where can semiconvection make a significant 
impact on the evolution of a star? 

As discussed in section semiconvective mixing only 
occurs in the presence of both stable composition gradi¬ 
ents and unstable thermal gradients, which are typically 
found in regions adjacent to convective zones. The most 
natural source of stabilizing composition gradients is nu¬ 
clear burning, and convective burning cores exist in main 
sequence stars (with M >1.2 Mq) as well as during core 
helium burning. 

In order to best understand the effects of our new semi¬ 
convective mixing prescription on stellar evolution, we 
first identify stars in which it may have a large effect. 
Figurej^shows a measure of the potential impact of semi¬ 
convection on the evolution of main sequence stars in the 
1 — 3 Mq range. We quantify this impact by compar¬ 
ing the predicted convective core sizes computed without 
semiconvection (and without smoothing, see above) un¬ 
der the Schwarzschild and Ledoux criteria, respectively. 
Because the Schwarzschild criterion does not take into 
account compositional gradients while the Ledoux crite¬ 
rion does, we expect the sizes of convective cores calcu¬ 
lated in the Schwarzschild case to be larger than in the 
Ledoux case in the presence of stabilizing compositional 
gradients. How much larger will depend on the size of 
outside the convective core, which in turn depends on the 
local nuclear burning rates and varies with stellar mass. 
Any semiconvective region, should it exist, must neces¬ 
sarily reside in between the Schwarzschild and Ledoux 
core boundaries. Its maximum extent is therefore well 
approximated by the difference between the two radii. 
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Fig. 3. — Measure of potential impact of semiconvection on the 
cores of main sequence stars of varying masses, as described in the 
main text. Core convection occurs for stars with M > 1.2 Mq, and 
the larger the computed ratio is, the larger the potential effect of 
semiconvection on the evolution of the star. A small window exists 
from M 1.2—1.7 Mq where there is a large enough compositional 
gradient outside the convection zone for semiconvection to have a 
significant impact on the star. See Figure]^ for an illustration of 
the compositional gradients in selected stars. 


which can be used as a proxy for the potential impact of 
semiconvection. 

Figure shows the time-averaged ratio of the differ¬ 
ence in convective core mass between the Schwarzschild 
and Ledoux cases, scaled to that of the Ledoux case, 
computed as 


AMcc \ 
Mcc,L / 



Mcc,L 


(26) 


where Mcc,s is the convective core mass obtained using 
the Schwarzschild prescription, Mcc,l is the convective 
core mass obtained using the Ledoux prescription, and 
the interval [rojTms] spans the period between the onset 
of core convection, and the end of the main sequence. 
Since this ratio is calculated from two different stellar 
models, the upper bound on the integral is taken to be 
the shortest lifetime - here the Ledoux model due to the 
smaller convective core. The advantage of this measure 
is that it does not depend at all on the semiconvective 
prescription employed, but merely probes which main se¬ 
quence stars have significant chemical gradients outside 
of their convective cores. Figure shows that the most 
favorable mass range to explore the effects of semicon¬ 
vection in main sequence stars is M « 1.2 — 1.7 Mq. We 
have performed these calculations on stars with masses 
up to 30 Mq without seeing another window where semi¬ 
convective mixing has a significant impact on convective 
core evolution during the main sequence. 

The significance of this particular mass range for semi¬ 
convection is best understood by looking at the composi¬ 
tional gradient generated by nuclear burning outside the 
convective core. Figure shows the difference between 
stars that are susceptible to semiconvection (1.3 Mq and 
1.5 Mq models), and a 2.8 Mq model which is not. The 
relatively weak temperature dependence of the proton- 


proton chain allows burning to occur in regions outside of 
the convective core at a rate which gradually drops with 
radius due to the temperature and density decreasing 
outwards. This radius-dependent burning rate causes the 
development of a composition gradient which steepens in¬ 
ward. While all stars in the mass range shown in Figure 
[^have nearly the same /i-profile outside their core for a 
given age, the position of the edge of the core is strongly 
dependent on the stellar mass and is located at radii 
with stronger or weaker /i-gradients. Lower-mass stars 
with smaller convective cores have larger ^-gradients just 
outside their cores, while larger-mass stars with larger 
cores have correspondingly weaker ^-gradients. This is 
why semiconvection can make the largest difference in 
the evolution with stars having the smallest convective 
cores. 

Stars may also be susceptible to semiconvection dur¬ 
ing core helium burning phases. However, helium burn¬ 
ing reactions are much more temperature sensitive and 
do not extend far out of the convective core. As a re¬ 
sult, there is not a large difference between the sizes of 
convective cores calculated under different convective cri¬ 
teria. A compositional gradient can also exist outside of 
the hydrogen burning shell due to the less temperature- 
sensitive proton-proton chain burning outside of the main 
CNO burning region, but this effect is much less dramatic 
than during the main sequence due to the higher temper¬ 
atures in the burning regions, and to the much shorter 
remaining stellar lifetime. We defer a discussion of the 
impact of semiconvective mixing on the later stages of 
high-mass stars to a future paper. 


4.2. Maintaining a composition gradient 

In order for semiconvective regions to persist for signif¬ 
icant periods of time, their turbulent mixing rates cannot 
be too large, otherwise the compositional gradient that 
causes them to exist will be destroyed. If this happens, 
the semiconvective region is simply converted into a con¬ 
vective region since it remains Schwarzschild-unstable. 
The efficiency of compositional mixing in semiconvective 
regions is determined by the size of the diffusion coeffi¬ 
cient Dsemi (see equation 25) and therefore by the layer 
height Hl - larger layer heights imply larger diffusion co¬ 
efficients. This suggests that stellar models with layered 
semiconvection can be split into three groups - (a) one 
in which the compositional mixing is too weak to have 
any influence on stellar evolution, (b) one in which the 
convective core size is modified, but the mixing is slow 
enough for the semiconvective regions to persist through¬ 
out the main sequence, and (c) one where the composi¬ 
tional mixing is fast enough to turn semiconvective zones 
into convective zones over a time scale much shorter than 
the main sequence lifetime. We anticipate the existence 
of a critical layer height, ML.crit) such that evolution with 
Ml ML,crit (case c) is similar to one without semicon¬ 
vection, but where the convective criterion is given by 
the Schwarzschild criterion. Meanwhile, if Ml ^ ML,crit 
(case a) we expect the star to evolve as if semiconvec¬ 
tion was absent, but the convection is determined by the 
Ledoux criterion. 

When semiconvection is enabled and material can be 
exchanged between the convective core and surrounding 
regions, there is a competition between nuclear burning 
which maintains the compositional gradient and semicon- 
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Fig. 4.— Profiles of composition gradients Vl (purple), radiative 
gradient Vrad (black), adiabatic gradient Vad (orange), and mass 
fractions of hydrogen blue), as a function of mass coordinate 

for stars that can have large semiconvective regions (1.3 Mq and 
1.5 Mq) as well as a 2.8 Mq star that cannot support significant 
semiconvection outside of its core. These profiles are obtained by 
evolving each star under the Ledoux criterion until the central mass 
fraction of hydrogen first drops below 0.4. The source of the mass- 
dependence in where semiconvection can occur comes from the size 
of the chemical gradient term in Vl immediately outside the con¬ 
vective cores. Large values of Vi, outside small convective cores in 
stars in the mass range 1.2 — 1.7 Mq cause a significant difference 
between the convective boundaries determined by the Ledoux and 
Schwarzschild criteria: the former is the position of the V l spike, 
located at m/MQ = 0.01,0.07, and 0.38, respectively. The latter 
is the position where Vrad = Vadi located at ralM(^ = 0.04,0.14, 
and 0.38, respectively. For higher mass stars such as the 2.8 Mq 
case shown, the compositional gradient outside the core is too small 
to allow for large semiconvective regions. 


vective mixing which tries to remove the gradient. We 
can therefore estimate i^L,crit by comparing the mixing 
timescale of the semiconvective region to the composi¬ 
tional evolution time scale in the convective core. The 
semiconvective mixing timescale is given by 

P . 

f . — _ semi ( 27 ) 

‘'semi / p. \ 5 J 

where /semi is the radial extent of the entire semicon¬ 
vective region and (Dsemi) is the mass-averaged diffusion 
coefficient in the semiconvective region. Similarly, we can 
compute the compositional change timescale as 


t 


/i,center 


/^center 

/^center 


(28) 


where /icenter is the mean molecular weight of the ma¬ 
terial at the center of the star. This is the same as the 
/r-evolution time scale for the entire convective core since 
it is fully mixed. Changing Hl affects the value of tsemij 
but not t/i,center- The Critical layer height is the value of 
Hl for which these two time scales are equal. Writing 
the diffusion coefficient as 



we can then estimate i?L,crit by equating tsemi and 


t 


/X,center ? 


H, 


L,crit 


H. 


D,. 


semi ,0 ^/x,center 


\ 2/3 


(30) 


In order to estimate /semi) we run a model under the 
Ledoux criterion and calculate where the boundary of 
the convective core would be under the Schwarzschild 
criterion; /semi is the distance between the actual Ledoux 
and hypothetical Schwarzschild convective boundaries 
Figurej^shows the resulting time-averaged iXn.crit values 
over the main sequence evolution of stars in the mass 
range of 1.2 — 3.0 Mq. Given that Hp « 10^° cm near 
the cores of such stars, we find that iXn.crit ranges from 
« 10^ - 10^ cm for stars in the mass range of greatest 
potential semiconvective impact, 1.2 — 1.7 Mq. 

While iXL.crit was determined from basic timescale ar¬ 
guments, we can easily compare it to the actual layer 
height s realized in numerical simulations. [Wood et al.| 
(20131 found that the minimum layer height is about 10 
times the scale of the basic instability, /qddC) with 


/oDDC = lOd = 10 


= 10 




.95 


d\nT 


dr 


dr 


.P9 


Pktv \ 

2<5 |V- VadlV 


1/4 


1/4 


10"* cm. 

(31) 


This gives a minimum layer height of « 10® cm, much 
larger than iXL,crit- This has a fundamental consequence: 
layered semiconvection is so efficient in main sequence 
stars that is accurately approximated by evolution un¬ 
der the Schwarzschild criterion, ignoring semiconvection 
altogether! 

Figure illustrates this statement with Kippenhahn 
diagrams of a 1.3 Mq star evolved with various mixing 
prescriptions. Although unphysically small, we can in¬ 
vestigate the effect of layered semiconvection using layer 
heights XXl ^ -flL.crit- As expected, we find that semi¬ 
convective regions are effectively radiative and the star 
evolves as if the Ledoux criterion was used. For the more 
realistic case where XXl ^ -ffL.crit, the evolution is nearly 
identical to the Schwarzschild case, as predicted above. 
Our results also highlight how sensitive the evolution is 
to the layer height. There are several orders of magnitude 
difference in Hp between models that are effectively the 
same as the Schwarzschild case and models that are effec¬ 
tively the same as the Ledoux case. Intermediate values 
of Hp initially have semiconvection zones of the same size 
that eventually transition into convective zones before 
core hydrogen depletion. Figure [^summarizes the con¬ 
vective core evolution for several different mixing criteria, 
showing convective cores that fall into three size groups. 
The smallest convective cores occur when smoothing is 
enabled, and semiconvective mixing is either turned off or 
is weak enough {Hp «C Lln.crit) that the evolution is effec¬ 
tively Ledoux. Convective cores are larger for the same 


® This is the same as the size of semiconvective regions in models 
where chemical transport is turned off (eg. by setting Dgemi = 0) 
or the layer heights used are small enough where mixing is incon¬ 
sequential. 
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Fig. 5. — Estimates of the time-averaged critical layer height, 
^L,crit 5 during main sequence evolution as a fraction of pressure 
scale height as a function of stellar mass. Semiconvective layer 
heights larger than -ffL,crit will rapidly mix the semiconvective 
regions before the star evolves off the main sequence, destroying 
the chemical gradient that supports them and turning the region 
convective. Smaller layer heights will allow layered semiconvec¬ 
tive regions to maintain their compositional gradient and survive 
through the main sequence. For all masses considered here, the 
critical layer height is orders of magnitude smaller than the ex¬ 
pected minimum layer height, which suggests that layered semi¬ 
convection cannot persist in main sequence stars and convective 
core evolution is therefore well approximated with that obtained 
under the Schwarzschild criterion. 

cases when smoothing is not used. Finaiiy, the iargest 
convective cores occur for efficient semiconvection (the 
oniy physicaiiy reaiizabie outcome, where ^ 7JL,crit) 
and is neariy identicai to the Schwarzschiid case. For 
these modeis, turning on smoothing does not have a no- 
ticeabie effect due to the strong mixing outside the con¬ 
vective boundary. 

Finaiiy, we show that the same resuits hoid in modeis 
with convective overshoot. The spikiness in the convec¬ 
tive core mass under Schwarzschiid evoiution can be re¬ 
moved with a smaii amount of convective overshoot. Fig¬ 
ure]^ shows the evoiution of a pure Schwarzschiid modei, 
a Scmwarzschiid modei with overshoot, and a semiconvec¬ 
tive modei with Hl ^ T^L.crit (thus effectiveiy as if they 
were run with the Schwarzschiid criterion). Semicon¬ 
vective modeis that are effectiveiy Schwarzschiid modeis 
remain so with convective overshoot, so the respective 
prescriptions combine without issue. 


5. CONCLUSION 

Convective core evoiution of main sequence stars in 
the mass range 1.2 — 1.7 Mq can be dramatically im¬ 
pacted by semiconvective mixing due the extended com¬ 
positional gradient outside the convective core coming 
from nuclear burning. We investigated the effects of the 
layere d semiconvection prescription given in [Wood et al. 
(20131 on the evolution of such stars, finding that there is 
a critical layer height, T^L.critj above which the evolution 
is effectively given by ignoring compositional gradients 
altogether and evolving the star using the Schwarzschiid 
criterion. For layer heights smaller than iLL.crit, the evo¬ 
lution is effectively given by ignoring the additional mix- 


Fig. 6.— Kippenhahn diagrams for the main sequence evolu¬ 
tion of a 1.3 Mq star using different mixing criteria. The pan¬ 
els on the left show the evolution of the convective core (blue) 
under the Schwarzschiid and Ledoux criteria. The convective 
core is much smaller in the Ledoux case because there is a sig¬ 
nificant stabilizing compositional gradient outside the core due 
to pp-chain burning. The main-sequence lifetime of this case is 
therefore significantly shorter. The panels on the right show the 
corresponding evolution with layered semiconvection using layer 
heights many orders of magnitude larger and smaller than the crit¬ 
ical layer height, ffn.crit ~ 10“® Hp. The core evolution of the 
larger Hp case is virtually the same as one obtained with a model 
using the Schwarzschiid criterion, while the evolution under the 
smaller Hp case is virtually the same as the Ledoux case. Only 
the Hp <C L^Ljcrit model shows large semiconvective zones (light 
blue) over the entirety of the main sequence, because the mixing 
is too slow to remove the compositional gradient. 


ing beyond the convective core and evolving the star us¬ 
ing the Ledoux criterion. This critical layer height is or¬ 
ders of magnitude smaller than the minimum layer height 
predicted from the underlying instability, so if layered 
semiconvection occurs within stars, we expect it to be 
very effective at mixing composition and quickly erasing 
the compositional gradient that allows it to exist. Such a 
star evolves as if the Schwarzschiid mixing criterion were 
employed. We also found that numerically smoothing 
the compositional gradient term can significantly change 
the sizes of convective cores. When using the Ledoux cri¬ 
terion, this V^-smoothing will artificially push the con¬ 
vective core boundary inwards, shrinking the convective 
core. This effect can be dramatic if no additional mixing 
is included beyond the convection zone, and can arti¬ 
ficially reduce the lifetimes of evolutionary phases with 
convective cores by up to fifty percent. 

It may be possible in principle to infer the sizes of con- 


of pulsation modes (Mazumdar et al.| 2006 Cunha & 

Metcalfe 

120071 ISilva Aguirre et al.^OlOaT 

IBrandao et al. 


mam sequence stars is difficult and the current KEPLER 
and CoRoT data for such stars is much more limited than 
for brighter objects such as red giants. Sample sizes of 
solar-like oscillators (both ma in sequence stars and sub 


Metcalfe et al. 

2014 

I, while those 

of red giants are in 

the tens of thousanc 

s ( 

Stello et al. 

2013 

|. The few in- 
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ferences of convective core sizes in main sequence stars 
show evidence of mixing beyond the Schwarzschild con¬ 
vection boundary, typically in terpreted as a constraint 
on the amount of overshooting ( Silva Aguirre et al.|2013 


Liu et al.|2014||Deheuvels|2015|). it therefore may not be 

possible to place constraints on the strength of semicon¬ 
vection in main sequence stars, so investigating its effect 
on evolved stars may be more fruitful. We are planning 
a future paper to examine this prospect. 

We thank Justin Brown, Chris Mankovich, and Bill 
Paxton for useful discussions. Code and inlists neces¬ 
sary to reproduce our models are hosted on mesastar.org. 
This work was supported under grants NSF AST 0847477 
and NSF AST 1211394. 


Fig. 7.— Main sequence evolution of the convective core mass, 
Mcc, for 1.5 Mq stars under various mixing schemes. Dashed lines 
indicate the same mixing parameters as solid lines, but with the 
default 7-point Gaussian V^-smoothing enabled. By adjusting the 
convection criterion, smoothing, and strength of layered semicon¬ 
vection, the convective core evolution falls into roughly three cases. 
The smallest convective cores occur when the Ledoux criterion (or 
layered semiconvection with with Hl ^fL,crit) is used in con¬ 
junction with V^-smoothing. Mid-size convective cores occur for 
the same cases but without V^-smoothing enabled. The largest 
convective cores occur with either the Schwarzschild criterion or 
with layered semiconvection where ^ -^^L.crit- For those cases, 
whether V^-smoothing is enabled or not does not have a significant 
effect since the mixing is strong enough to destroy compositional 
gradients. 
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Fig. 8.— Kippenhahn diagrams for the main sequence evolution 
of a 1.4 Mq star showing the effects of overshoot on smoothing out 
the temporal evolution of the convective core. Dark blue regions 
are convective, light blue are semiconvective, and orange corre¬ 
spond to overshooting. The top panel shows the evolution under 
the Schwarzschild criterion. The middle panel shows the same evo¬ 
lution but also with overshoot turned on (using / = /q = lO""^). 
The bottom panel shows evolution with layered semiconvection at 
a value of = 10~® Hp and the same overshooting parameters. 
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